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ABSTRACT 

We investigate the evolutional and environmental effects on star formation efficiency 
for more than 400 merging galaxies. The ~400 merging systems, with photometric red- 
shifts smaller than 0.7, are obtained from a catalog of ~15000 morphologically identified 
merging galaxies derived from observations of the Canada-France-Hawaii Telescope. We 
also obtained the IR data of the merging galaxies from the Spitzer Wide-area InfraRed 
Extragalactic Survey (SWIRE). The redshift differences Az between the member galax- 
ies of these merging pairs show a large distribution with < Az < 0.4. We divide our 
merging pairs into two sub-samples with Az < 0.05 and > 0.05 for further analyses. 
We find a statistically significant anti-correlation between the specific star formation 
rate (SSFR) and the separation of the merging galaxies for both sub-samples. Our 
analyses also show that although most of the merging systems do have enhanced star 
formation activity, only very rare ones display extremely high SFRs. Additionally, the 
SSFR of the merging galaxies also decreases when the magnitude difference between 
two member galaxies becomes large. However, we find that for the merging pairs with 
large luminosity contrast, the fainter components show higher SSFR than the brighter 
ones. Finally, there is a higher fraction of gas-poor mergers in galaxy clusters, and the 
SSFR of gas-rich mergers is reduced in cluster environments. 

Subject headings: galaxies: general — galaxy: interactions — infrared: galaxy 



1. Introduction 

Galaxy interaction and merging play an essential role in galaxy formation and evolution. Nu- 
merical simulations sho wed that merg i ng pr ocesses might transform blue late-type galaxies into 
red early- type galaxies (jTissera et al.l l2002l ). It was also predicted that tidal interactions be- 
tween galaxy pairs could transform spiral and irregular galaxies into SO arid elli ptical galaxies 



(Toomre &: Toomre 



1972 



Mihos Hernquist 



1996 



Barnes &: Hernquist et al.lll996l ). S ome obser- 



vatio i is showed that ellipticals forni ed through the merger of spirals (e.g., the Antennae: ISchweizer 



1982 



Hibbard van Gorkoml 119961 ). Non-axisymmetric gravitational potential arising during the 



galaxy interactions might produce tidal torques that can induce strong gas infiows, which might 
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trigger an intense burst of star formation. Galaxy mergers are therefore w e ll known to be re- 
sponsible for trigge ring enhanced star formation activity (e.g.. lKennicuttlll998l : iBarton et al.ll200d : 
Lambas et al.ll2003l ). 



An interesting issue of merging galaxies is how the triggered star formation activity is related 
to the internal physical parameters of the galaxy mergers, such as pair separations and mass ratios. 
It has been reported that star form ation rate (SFR) i s anti-correlated with the pair projected 
separation and radial velocity (e.g., iBarton et al.l 120001 : iLambas et al.l l2003l ) . Since SFR scales 
with galaxy mass, it is better to use the specific star formation rate (SSFR), defined as the SFR 



divided by the stell ar mass, to denote star formation activity of galaxy (e.g.. iNikolic et al.ll2004l ). 
Nikolic et al.l (|2004l ) discovered that the SSFR also shows an anti-correlation with the projected 
separation. In addition, the induced star formation activity is also related to the mass ratio of 
the galaxy pair; the pairs with commensur ate luminosities trigger stronger star for mation activity 
than tho se with large l u minos ity contrast (jLambas et al.l l2003l : IWoods et al.l |2006| ) . In the minor 
mergers, 



Lambas et al. (j2003l ) found that brighter membe rs show a highe r prob ability to have 



tidally enhanced star formation than the fainter ones, but Woods &: Geller (j2007l ) reported that 
the fainter galaxies show higher SSFRs than the brighter ones. 

Besides th e internal par a meter s, the triggered star formation activity might also be affected by 
environments. lAlonso et al.l (|2004l ) found that the pairs in clusters required a smaller separation 



than those in the fields to show obvious enhanced star formation activity by studying 2-degree field 
(2dF) pairs. They also pointed out that the pairs in dense environments showed lower enhanced 
star formation activity than those in the fields. Nevertheless, it is difficult to determine whether 
this SFR inhib ition is true or mainly results from the effect of the morphology-density relation 
(jDresslerlllOSd ) since ellipticals, which are dominant in clusters, have much smaller star formation 
activity than spirals. 

Early studies based on UV/optical observations were susceptible to dust extinction (e.g. 



Barton et al 



200C: 



Woods fc Cxellei 



Lambas et al 



200 



I" 



2003 



Nikolic et al 



2004 



Alonso et al 



2004 



Woods et al. 



2006 



During the merging process the gas infiows induced by tidal torques of 
galaxy interaction are converted into new stars with large amounts of dust, which is heated by 
the young stars and radiates energy at mid-infrared (MIR) to far-infrared (FIR) wavelengths; this 
indicates that infrared (IR) is an ideal tool to study mergers. 

Furthermore, previous studies on the relevant topics were all based on merging galaxies selected 
using the inter-galaxy separations. Such methods would tend to select interacting galaxies at 
relatively early stages of the merging processes. To statistically study the influence of the merging 
process and environments, we need a large un-biased merging galaxy sample based on consistent 
searching criteria. In t his study, we select a sample of morphologically-identified merging galaxies 
dHwang fc ChaneJbood l. which are not prone to specific merging stages, and we combine the infrared 
and optical data to investigate the properties of these merging galaxies. 



This paper is organized as follows. Section 2 introduces our data and data processing. In 
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Section 3, we discuss the photometric redshifts of the galaxy pairs and estimate the star formation 
rates. In Section 4, we study the distribution of SSFR in the merging galaxies. The relationships 
of the SSFR with the projected separation and the mass ratio of the galaxy pairs are analyzed in 
Section 5. In Section 6 we discuss the environmental effects on the SSFR of the mergers. Finally, 
Section 7 summarizes our main findings. Where needed, we assume a standard ACDM cosmology 
with Q,M = 0.3, Q,A = 0.7, and Hq = 70 km s^^ Mpc^^ throughout this study. 



Data and Data Reduction 



2.1. Sample of Merging Galaxies 



The merging galaxies in this study are obtained from a catalog of interacting and merging 
galaxies derived fro m the Red Sequence C luster Survey 2 (RCS2) of Canada-France-Hawaii Tele- 
scope observations (jHwang &: ChaneJ l2009l ) . This catalog is the largest catalog of galaxy-merging 
systems assembled to date and includes more than 15,000 interacting and merging galaxies via 
morphological classification. Possible candidates were first selected by an automatic morphological- 
pattern-recognition package and subsequently inspected by human eyes to assure their credibility 
in the final catalog. These sources constitute a reliable and uniform sample of merging galaxies 
for further study. In addition, the RCS2 data used to develop this catalog covers extensive and 
separate areas up to 422 deg^, providing a great opportunity to conduct multi- waveband researches 
combined with o ther large surveys, su ch as the Spitzer Wide-area Infrared Extragala ctic Legacy 



Surv ey fSWIRE: iLonsdale et al.ll200l 2004) and the Sloan Digital Sky Survey fSDSS: lYork et al 
2OO0l ). 



The SWIRE is the largest extragalactic Spitzer Legacy Science programs, imaging ~49 deg^ 
at the wavelengths of Infrared Array Camera (IRAC) 3.6, 4.5, 5.6, 8 /xm and Multiband Imaging 
Photometer (MIPS) 24, 70, 160 nm. The SWIRE project has mapped seven high-latitude fields; 
the RCS2 regions are overlapped with 3 of these 7 fields — ELAIS Nl, ELAIS N2, and Lockman 
Hole. These 3 SWIRE fields are listed in Table 1. We note that these fields have also been 
specifically selected to be far away from the ecliptic to avoid possible contaminations from wandering 
asteroids as well. With in these 3 fields there are 540 merging systems identified in the catalog of 
Hwang fc ChaMi tood ). 



2.2. Optical Data from SDSS 

In order to have accurate photometric information for these morphologically identified merging 
galaxies, we also obtain the optical data from the SDSS survey. The SDSS is a combined imaging 
and spectroscopic survey of 10^ deg^ in the north Galactic cap and a smaller region in the south. 
The SDSS imaging gathers data in five broad bands, ugriz, ranging from 3000 to 10,000 A. We use 
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the imaging data in the seventh SDSS data release (DR7) in this study. 

We first retrieve a hst of galaxies within an angular distance of 10" from the 540 merging 
systems. We visually examine each image of the 540 merging systems to identify their SDSS 
counterparts. It turns out that a small part of the merging systems are identified as foreground 
stars by SDSS. We find 60 merging systems composed of at least one component identified as stars 
by SDSS; we reject these sources from our sample. We then classify the remaining 480 merging 
systems, totaling to 723 galaxies, into 3 groups according to the number of their members: merger 
(MG; 1), close pair (CP; 2), and close multiple (CM; more than three). The final sample consists of 
253 MGs, 213 CPs, and 14 CMs. In the following analyses, we only focus on the MGs and CPs to 
study the relations between the SSFR and some physical parameters, such as pair separation, pair 
magnitude difference, and the environments. The CMs are ignored in these analyses because they 
are complicated and too difficult to quantify the effects of the physical parameters. Some examples 
of these merging systems are shown in Figure 1. 

The SDSS database contains various measured magnitudes for each detected source; we se- 
lect model magnitudes for our targets. We also acquire the photometric redshifts for our galaxies 
because of their lack of spectroscopic observations. Notwithstanding the accuracy of the photo- 
metric redshifts is not as good as that of spectroscopic redshifts, the former still provides sufficient 
information for extragalactic and cosmological researches, especially in the statistical study of the 
evolutionary properties for faint objects. 



2.3. IR Data from SWIRE 



We search the IR counterparts of the merging galaxies in the SWIRE fields to study their 
properties of star formation activity. We mainly utilize the SWIRE data at 3.6 and 24 fj,m. The 
3.6 /im luminosity can be used as a tracer of the stellar component since the NIR emission is 
dominated by o l der stellar populatio n s and less susceptible to recently formed massive stars (e.g., 



Wu et al. 



2005 



Davoodi et al 



20061 : 



Li et al. 2007: Hancock et al 



20071 : 



Rowan-Robinson et al 



20081 ). On the othe r hand, the MIR emission can be used t o estimate the total IR l uminosity for a 



varie ty of galaxies (IChary k Elbadl200ll : lElbaz et al.ll2002l : IPapovich k Belli I2OO2I : iTakeuchi et al 



2OO5I ). We use the 24 nm flux to derive the total IR luminosities and thus the SFRs of the merging 
galaxies. 



2.3.1. IR Counterparts from SWIRE Catalogs 

We obtain the IR data for bright sources from the SWIRE catalogs. We first check the SWIRE 
flag — extended source flag — for each source in SWIRE catalogs to exclude "definitely point-like" 
objects. The processed SWIRE catalogs are subsequently cross-correlated with our 723 merging 
galaxies; the IR counterparts are determined by identifying the closest SWIRE sources within an 
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angular distance of 5" from each merging galaxy. We find 482 and 152 galaxies with IR counterparts 
at 3.6 and 24 //m respectively. The SWI RE catalogs provide several kinds of flux measurements; 
for our objects, we select the Kron flux (jKronl 119801 ) . which is especially preferable for extended 



sources. 



2.3.2. IR Counterparts from SWIRE Images 

The SWIRE catalogs only contain sources with sufficient fluxes. For example, the threshold 
for the detection at 3.6 /im required SNR > 10, which is equivalent to 10 //Jy for almost all of the 
survey area. To retrieve faint sources which were not included in the SWIRE catalogs, we extract 
the fluxes within apertures of 1.9" (~1.6 FWHM) and 5.25" (~1.9 FWHM) centered at the position 
of each merging galaxy in the IRAC-3.6 fim. and MIPS-24 ^m images separately. These images 
are downloaded from the NASA/IPAC Infrared Science Archive and have been post-processed and 
co-added. They are the same ones used to make the SWIRE catalogs. We obtain the IR fluxes 
of the merging galaxies if the corresponding extracted objects have an SNR greater than 3 sigma. 
Combining the faint sources and the objects from the SWIRE catalogs, we find 600 and 621 out 
of 723 merging galaxies with IR counterparts at 3.6 and at 24 ^m respectively. There are fewer 
galaxies with IR counterparts at 3.6 than at 24 /xm because some of the galaxies are out of the 
IRAC imaging fields. 

Since we can not detect very faint sources at high redshifts, we might have some biases in 
selecting the merging pairs. For examples, if we divide our sample into major merging pairs with 
Am^ < 1.5 and minor merging pairs with Am^ > 1.5, we might tend to detect the minor mergers 
at relatively low redshifts becasue we can have more detectable low-luminosity sources at lower 
redshifts. Nevertheless, the redshift distributions of the major and minor mergers are not distinct; 
the mean redshift for the major mergers is 0.17 it 0.11 and that for the minor mergers is 0.15 it 0.10. 
This indicates that this flux-limited bias is negligible. 



3. Sample Characteristics 

3.1. Photometric Redshifts of Close Pairs 

The redshift distribution of our sample is shown in Figure 2. The farthest merging galaxies 
could reach out to redshifts as far as z ~ 0.7; the observed 24 /im flux for this largest redshift corre- 
sponds to the rest-frame 14.1 /im flux, still within the MIR range. Figure 3 shows the distribution 
of the redshift difference, Az, between two member galaxies for each CP. At first glance, this wide 
distribution might imply that our sample is contaminated by spurious galaxy pairs; i.e., some pair 
galaxies might not be really close in physical space. However, we note that our sample is selected 
through morphological identification; most CPs exhibit clear features of merging galaxies, such as 
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"tail" and "bridge" . Some examples of the CPs with Az > 0.05 are shown in Figure 4. As can be 
seen, these CPs do show the features of merging galaxies in spite of their large redshift difference. 
We thus conclude that the redshift differences for most of these CPs are not caused by projection. 

The discrepancy of photometric redshifts in some of the CPs is not surprising. Enhanced 
star formation in the interacting galaxies may affect their color caused by strong line emission 
and dust absorption in the star forming region. This would influence their photometric redshift 
estimates derived by fitting observational data with synthetic templates. We have considered 
different criteria in obtaining a better redshift estimate between two interacting galaxies, such as 
selecting the photometric redshift from the brighter galaxy or selecting the smaller redshift. We 
find that our derived results are all statistically consistent in spite of the selecting criteria used; 
one of the main reasons is that the redshift effects are substantially canceled out when we consider 
the SSFRs derived from the ratios of luminosities (fluxes) of different wavebands. For simplicity, 
we assume the smaller photometric redshift between two member galaxies as the true redshift for 
those CPs with > 0.05 in presenting our results. 



3.2. Estimating Total IR Luminosity and Stellar Mass 



Emission of MIR is an ideal tracer of the total IR luminosity. Previous studies have developed 
libraries of luminosity- o r color-dependent galaxy templates to calculate the total IR luminosity 



from 24 flux densities (jPale et al 



2001 



Charv Elba2 



2001 



Dale Helou 



2002 



Lagache et al 



20031) . We derive the t otal IR luminosities from the 24 //m flux densities adopting the templates from 



Chary &: Elbaa (120011 ): these templates consist of a combination of 105 spectral energy distributions 



(SEDs) from normal to starburst galaxies. 
We employ the SDSS z-band (Ag// 



3.6 fim flux density was usually used as an estimate for the stellar mass (e.g. 


Wu et al. 


2005; 


Davoodi et al. 


2006; 


Li et al. 


2007; 


Hancock et al.l 2007: 


Rowan-Robinson et al. 


2008). As found 



bv ICharlotI (jl996l ) and iMadau et al.l (]l998l ). the mass-to-NIR ratio, independent of either galaxy 
color or Hubble type, is less susceptible to the star formation history and the same for all types of 
galaxies. Nonetheless, there are 123 sources in our sample without information at 3.6 /im due to 
detection limit or to the lack of IRAC observations. On the other hand, it ha s also been reporte d 
that z-band luminosity can be used as an un-biased estimate of the stellar mass (iNikolic et al.ll2004l ) . 
To check the credibility of using the SDSS z-band as the stellar mass indicator in our sample, we 
plot the 3.6 fxm versus the SDSS z-band luminosity for the 600 sources with both 3.6 ^um and 
z-band detections in Figure 5. The 3.6 urn and the S DSS z-band luminosities are calibrated to the 
rest frame using the templates from lCharv &: Elbaa (|200ll ). As can be seen, the 3.6 ^m luminosity 
is proportional to the SDSS z-band luminosity. The z-band luminosity provides a reliable estimate 
of stellar mass. We have hence obtained an estimate of the stellar mass based on the SDSS z-band 
luminosity for every galaxy in our sample. 
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Throughout this study, we use Lm/Lz to represent the SSFR. The SFR scales with the galaxy 
mass and might fail to reflect the star formation efficiency; in contrast, the SSFR, which has taken 
into account the effect of stellar mass, is a more suitable estimate of star formation activity in 
galaxies. The LjrIM^ ratio has been employed to represent the SSFR in previous studies (e.g., 



Lin et al.l 120071 ) . They used rest-frame (B-V) colors and absolute S-band magnitude to derive 



stellar mass. We note that th e .B-band i nagnit ude is susceptible to dust extinction and requires 
careful calibration as done by iLin et al.l (|2007l ). In this paper, we use rest-frame SDSS z-band 
luminosity to represent the stellar mass because the SDSS z-band luminosity is proportional to the 
rest-frame IRAC 3.6 //m luminosity (Figure 5), which is less susceptible to star formation history 
and dust extinction and is a better mass tracer. For each CP, we use Lm/Lz to represent the total 
SSFR by summing the total IR luminosity of both galaxies over the total z-band luminosity of both 
members. We note that the distance influence on the SSFR is small because the division of Lm 
by Lz will have eliminated most of the distance effects. Thus, the uncertainties in the photometric 
redshift estimates are not as vital as might have been expected. 



4. SSFR Distribution of Galaxy Pairs 

We first explore characteristic features of the SSFR distribution of the CPs. Upper panel of 
Figure 6 shows the SSFR distribution of 213 CPs. This distribution shows a significant dip at 
Ljji/Lz ~ 1 — 1.5, which looks like a boundary of two diffe rent distributions. T o understand the 



origin of this boundary, we adopt the SED templates from iPolletta et al.l (j2007l ) to generate the 
Ljr/Lz values for various types of galaxies. Table 2 summarizes the values of Ljr/Lz for a variety 
of galaxies ranging from ellipticals to ultra-luminous infrared galaxies (ULIRGs). As can be seen 
in Table 2, while the galaxies with Lm/Lz below unity are normal ones, the galaxies with Lm/Lz 
larger than unity all exhibit enhanced star formation activity. This indicates that the bimodal 
distribution of Figure 6 actually symbolizes the galaxies with and without enhanced star formation 
activity. Some galaxies of the CPs are influenced by tidal interactions and manifest enhanced star 
formation activity whereas others still remain unaffected. 

We try to find suitable distribution functions to describe the SSFR distribution. First, the CPs 
on the left side looks like a normal distribution, which seems to be reasonable for normal galaxies. 
The probability density function (p.d.f.) of a normal distribution is 

1 -{x-lJ.f 

fix\ fi,(7) = ^=e , (1) 



where and a are the mean and standard deviation of the variable x. On the other hand, we 
might expect the distribution of galaxies with enhanced star formation activities to be lognormal; 
that is, only very few sources have extremely star formation activities. The right hand side of the 
SSFR distribution shows a concentration of the Ljr/Lz values at 1.5-2 and a long tail toward the 
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upper end, also corresponding with the characteristics of a lognormal distribution. The p.d.f. of a 
lognormal distribution is 

1 —tlnx — y,)'^ 

f{x\^^,a) = e^^, (2) 



where fi and a are the mean and standard deviation of the variable's natural logarithm. The 
lognormal distribution is closely related to normal distribution. If x is distributed lognormally with 
parameters // and a, log{x) is distributed normally with the mean fi and standard deviation a. We 
adopt a normal function to fit the left part and a lognormal distribution function to fit the right 
part of the SSFR distribution. Table 3 lists the derived parameters, and the results are plotted in 
the middle panel of Figure 6. 

The lognormal fits are somehow unsatisfactory. This is due to the fact that at the low end there 
are always galaxies with "normal" star formation efhciency, and on the high end, there arc more 
galaxies than model predicted. In order to correct these shortcomings, we construct a modified 
lognormal p.d.f. : 



f{x \^x,a) = {x- If e-^^ . (3) 

xav z-K 



We note that such modification is purely from mathematical reasoning. The derived parameters 
are also listed in Table 3, and the fitting results are plotted in Figure 6 (lower panel). The data 
are better fitted with a mixture of the normal and the modified lognormal distribution functions. 
We note that the modified lognormal distribution is still very close to the standard lognormal 
distribution. We conclude that the SSFRs of the CPs that remain unaffected by tidal interactions 
are normally distributed, but the SSFRs of those with enhanced star formation will approximately 
follow a lognormal distribution. 



5. Dependence of SSFR on Physical Parameters 

5.1. Projected Separation 

We explore effects of the projected separation Ar^ of the merging galaxies on the SSFR. We 
ignore the CPs without 24 /im detection; the total remained CPs in this analysis is 181. Figure 7 
shows the dependence of L/^j/L^ on the projected separation; the circle and plus signs represent the 
CPs with < 0.05 and Az > 0.05 respectively, and the triangle represents the MGs. In Figure 
7, the Ljr/Lz shows a generally declining envelop as a function of Ar^. To test the statistical 
significance of this relation, we compute their correlation coefficients and probabilities. The linear 
correlation coefficient r for the overall sample is —0.18, with the probability of the correlation 
P = 0.98, indicating Lir/Lz and Ar^ are significantly anti-correlated. The sample with Az < 0.05 
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shows a correlation of r = —0.27 with P = 0.94, and the sample with Az > 0.05 shows a correlation 
of r = —0.16 with P = 0.94. To reduce the effect of outliers in the correlation test, we also apply 
the Spearman rank method, which is more robust than the linear correlation. A Spearman rank 
correlation test also confirms the anti-correlation between the SSFR and the projected separation; 
we obtain = —0.17 and Pg = 0.97 for the overall sample, = —0.33 and Pg = 0.98 for the 
sample with Az < 0.05, and = —0.17 and Pg = 0.95 for the sample with Az > 0.05. All these 
tests show that the anti-correlation between the SSFR and the observed projected separation Avp 
of merging galaxies is statistically significant. 

Our results qualitatively agree with previous studies. It was reported that SSFR is anti- 



corre lat ed with pa i r pro jected separation for the merging galaxies of local universe (jNikolic et al 



20041 ). iLin et al.l (j2007l ) also found a declining envelop of L/r/M^, as a function of projected 
separation for pairs at higher redshift (z ~ 0.1-1.1); however, no statistical analysis was provided 
and is difficult to evaluate and to compare with their results quantitatively. 

In our sample, the CPs with similar projected separations spread a wide range of the Lm/Lz 
values. We note that CPs with similar projected separations might in fact have very diverse real 
separations; this could cause large scattering in the Lm/Lz values. CPs with similar physical 
separations but undergoing different merging stages would also contribute to the scatter of SSFR 
because some CPs may just start to approach, while others might have experienced merging pro- 
cesses more than once and had drained most of the gas in the early processes. Likewise, the same 
fact might explain why the L^/Lz values of MGs also scatter over a wide range. It was also 
noted that the SSFR is related to the g alaxy types or intrinsic properties of the merging galaxies 



(Lambas et al. 



2003 



Nikolic et al.ll2004l ) and the luminosity contrast (or mass contrast) between 



member galaxies. Given all these uncertainties, our results still show a significant anti-correlation 
between the SSFR and the observed projected separation Ar^ of merging pairs; this demonstrates 
that the internal separations between merging galaxies must have strong influence on the SSFRs. 



5.2. Mass Ratio 

We investigate the dependence of the SSFR on the mass ratio of our CPs. We use the z-band 
magnitude difference between the pair galaxies, Am^, to represent their mass ratio. We divide our 
sample into major merging pairs with Am^ < 1.5 and minor merging pairs with Am^ > 1.5. The 
number of the major merging pairs is 134, including 38 CPs with Az < 0.05 and 96 CPs with 
Az > 0.05; the 47 minor merging pairs contain 11 with Az < 0.05 and 36 with Az > 0.05. The 
distribution of the z-band magnitude differences for the CPs is shown in Figure 8. 

We first examine the relation between Lm/Lz and Avp to test the influence of the luminosity 
contrast for the major and minor merging pairs. Figure 9 shows Lm/Lz versus Arp of these two 
types of merging pairs. We find a significant anti-correlation between the Lm/Lz and Arp for the 
major merging pairs with = —0.20 and Pg = 0.98 using the Spearman rank test. We further 
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divide the major merging pairs into two subgroups with Az < 0.05 and Az > 0.05 separately. The 
statistical significance of the anti-correlations of these two subgroup are = —0.43 and Ps = 0.99 
for the major merging pairs with Az < 0.05 and = —0.18 and Ps = 0.92 for the ones with 
Az > 0.05 respectively. In other words, the correlation between Ljji/Lz and Arp becomes weaker 
for the major merging pairs with Az > 0.05; this might be caused by the errors of the true distance 
determination for some sources in this sample. 

In contrast, the minor merging pairs display no correlations between Lj^/L^ and Arp. The 
Spearman test show = —0.08 and Pg = 0.41 for the minor merging pairs. There are no anti- 
correlations f or both subgroups w ith Az < 0.05 and Az > 0.05 either. This is consistent with 
the results of IWoods et al.l (|2006l ) and implies that the intrinsic properties of the minor merging 
galaxies might be dominant factors in determining the SSFR of the minor merging galaxies. 

We next study the relation between the SSFR and the magnitude difference in these merging 
galaxies pairs. Figure 10 shows the relation between Ljr/Lz and Auiz for our sources. The 
results of the Spearman rank test exhibit a clear anti-correlation between the Ljji/Lz and Am^ 
with = —0.21 and Ps = 0.99 for the overall sample; the correlations are slightly weaker with 
= —0.30 and Ps = 0.97 for the pairs of Az < 0.05 and = 0.16 and Pg = 0.93 for the pairs of 
A z > 0.05. Th i s anti -correlation between the Ljji/Lz and Am^ is also consistent with the result 
of lWoods et"ZI (|200d ). 



We also compare the SSFRs of the faint galaxies and the bright galaxies in the minor merging 
pairs. To estimate the stellar mass of galaxies, we apply the relation between the z-band and 
3.6 fim luminosities of our sample as shown in Figure 5, 



log 







(0.17 ± 0.13) + (0.92 ± 0.01)/off ( ^ 



and the relation found bv lLi et al.l (|2007l ) 



log[:-:^]= (1.34 ± 0.09) + (1.00 ± Qm)log ( ^-^'-[3.6^m] \ 



Mr. 



(4) 



(5) 



We note that equation (4) also holds in different redshift ranges, indicating that its redshift depen- 
dence is small. 

Dwarf galaxies are excluded in calculating the SSFRs since the dwarves might have depleted 
their gas and their SSFRs are not affected by tidal interaction anymore. We exclude dwarf galaxies 
with log(M/MQ) < 9, which corresponds to galaxies about o ne order magnitude fainter than the 
characteristic luminosity of field galaxies ()Blanton et al.ll2005l ). 



We calculate the individual SSFRs of the faint galaxies and the bright galaxies in the minor 
merging pairs. As shown in Figure 11, all the faint galaxies have higher SSFRs than the bright 
galaxies. Their average SSFRs are listed in Table 4; the average SSFRs of the major merging 
galaxies are also listed as a comparison. Although the major merging pairs have higher SSFRs 
than the minor pairs, the average SSFR of the faint galaxies in the minor merging pairs is much 



- 11 - 



higher than that of the major merging pairs. This is due to the fact that the SSFRs of the bright 
galaxies in minor merging pairs are much lower than those of the major merging pairs; the very low 
SSFRs and the high mass of the bright galaxies compensate the high SSFRs of the faint galaxies in 
the minor merging pairs. As a result, the major merging pairs have higher SSFRs than the minor 
ones on average. Our results indicate that the SSFR of faint galaxies are more susceptible to the 
effects of tidal interaction. 



Th is result is in agreem e nt with that o f IWoods &: Gelleri (|2007l ) but are inconsistent with 



that of Lambas et al. ( 2003 ). Lambas et al. ( 20031 ) found that in the minor merging pairs, the 



effects on star formation activities are more important in the brighter galaxies. We note that their 
faint galaxy sample including sources down to log(M/MQ) ~ 8, which might contain many dwarf 
galaxies. These dwarf galaxies might have depleted most of their gas and could not trigger any new 
star formation even with very strong tidal interaction. This shows that the intrinsic properties of 
galaxies have to be taken into account when we consider the effect of tidal interaction in a merging 
pair. 



5.3. Estimation of Superposition 

The redshift differences Az between the member galaxies of our merging pairs show a large 
distribution with < Az < 0.4. It is difficult to directly quantify the contamination because 
our sources have been specifically selected via morphological pattern recognition. This selection 
criterion has largely removed superposed false pairs. Nonetheless, we might evaluate the effects of 
overlapping pairs on our results from some correlation tests. We note th at the internal sepa rations 



between merging galaxies have a strong correlation with the SSFRs (e.g.. lNikolic et al.ll2004l ). so the 
correlation results could reflect the reliability of our sample. We divide the pairs with Az > 0.05 into 
three subgroups according to their redshift differences and calculate the Spearman rank correlation 
test on the SSFR and the separation for each subgroup, as listed in Table 5. 

We find that there are no correlations between Lm/Lz and Arp for CPs with 0.05 < Az < 0.2 
but a significant one for CPs with Az > 0.2 (r^ = —0.24 and Pg = 0.98). In Figure 3, we find a small 
bump around 0.1 < Az < 0.2 and the distribution becomes flat for Az > 0.2. From the correlation 
tests, we speculate that the bump around Az ~ 0.1 to 0.2 might be caused by superposition, whereas 
most of the sources with Az > 0.2 are actually true pairs with wrong redshift determination. This 
unusual result is caused by the fact that there are only a few high-redshift sources and thus the 
probability of superposition involving one member at high redshifts is low. It turns out that there 
are very few overlapping pairs with true large Az. From the bump of Figure 3, we estimate that 
there are about fifty contaminating sources by comparing the distribution of large Az sources. 
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6. ENVIRONMENTAL EFFECTS ON SSFR 

6.1. Environmental Distribution of Wet /Dry /Mixed Merging Pairs 

Environments also have strong influence on the merging galaxies. To study the environmental 
effects on galaxy evolution, we compare the distribution of merging galaxy types in the fields with 
that in the clusters. The c luster regions for comparison are selected from the maxBCG cluster 



catalog (jKoester et al.ll2007l ). which covers similar regions of our sample galaxies and is the largest 
sample of observed galaxy clusters available by far. The criteria of assigning a CP in a cluster are 
(1) the redshift difference between the CP and its host cluster is smaller than 0.05, and (2) the 
projected separation between the CP and the center of the cluster is smaller than 1000 kpc. 

We classify the types of galaxies according to their SSFRs. A galaxy is classified as a red (or 
dry) galaxy if Lju/Lz < 0.5 and a blue (or wet) one if Lf^/Lz > 0.5. Since some galaxies have 
no 24 fim detection with 3a upper limit, the galaxies with the Lm/Lz upper limit lower than 0.5 
are also classified as red galaxies. There are about 16 sources, which have the Ljr/Lz upper limits 
larger than 0.5 and are difficult to classify. We could either assign these uncertain sources as blue 
galaxies or red galaxies. In the former case, we would have the upper limit for the total number 
of the blue galaxies and in the latter case the upper limit for the red ones. The results from the 
former and latter classifications are represented as class 1 and class 2 in Table 7 separately. 

As can be seen in Table 7, the dry-to-wet ratio in the cluster is higher than that in the field. 
The dry-to-wet ratio of class 1 serves as a lower limit of the ratio because those galaxies with 
the LiFi/Lz upper limit higher than 0.5 might actually be red galaxies. On the other hand, the 
dry-to- wet ratio in class 2 is the upper limit one as we classify all galaxies with the Liji/Lz upper 
limit into red galaxies. Hence, the dry-to-wet ratio in the clusters is around 0.25-0.75 and the ratio 
in the field is around 0.03-0.09. In addition, only about 3% of all the wet merging galaxies are 
located in clusters, but about 24-29% of dry merging pairs and 18-21% of mixed merging pairs 
are in clusters. Dry mergers show a higher tendency to populate the clusters than other types of 
mergers; this tendency might be related to the morphology-density relation in which early-type 
galaxies preferentially populate h igh density eri vironments while late-type galaxies are likely to 



inhabit low density environments (|Dressleijll98d ). 



6.2. Inhibition of Star Formation Activity in the Cluster 



The infiuence of cluster environment on the SSFRs of galaxies is obvious for our merger sources 
(MGs). These interacting sources have merged together and the individual galaxies can not be 
identified. We find that the average SSFR of the MGs in the fields (156 sources) is 2.47 ± 0.20, 
and the average SSFR of those in the clusters (9 sources) is 1.86 = fc 0.41. We not e that, however, 
we could not exclude the effects of the morphology-density relation (|Dresslerill980l ) on the different 
SSFRs of these MGs since we do not know the galaxy types of the progenitors of these mergers. 
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In order not to confuse the environmental effects with intrinsic pair properties, we further 
select only the wet major merging pairs both in the clusters and in the fields to explore whether 
their SSFRs are different. We use Ljr/Lz = 0.5 as the dividing point to separate the wet galaxies 
from the dry ones. The result is shown in Table 8. 

We find that the SSFRs of the galaxy pairs are reduced in the clusters. The average SSFR of 
the wet major merging galaxies in the clusters is 2.28 and that in the fields is 4.36. We note that 
the average separation of the pairs in the clusters is 13.29 kpc, similar to that of pairs in the field, 
13.14 kpc; this suggests that the separation effect is not important. Besides, since we consider only 
wet major me rging galaxies , the d ifferent SSFRs can not be associated with the morphology-density 
relation as in lAlonso et al.l (|2004l ). 



Furthermore, we find an anti-correlation between the SSFR and Lz among the MGs in the 
fields (Figure 12). The Spearman rank test shows a highly significant result with = —0.23 and 
Pg = 0.99, indicating that in the fields the enhanced star formation activity of the more massive 
MGs is weaker than that of the less massive ones. Nonetheless, there is no correlation for the MGs 
in the clusters with Pg < 0.50. This implies that the star formation activity and/or merging history 
are very different for the MGs in the fields and those in the clusters. 

Our results demonstrate that the environments of clusters do inhibit the star formation activity 
among galaxy pairs. This inhibition might be due to (1) ram pressure stripping of c old interstellar 



medium of galaxies falling into the cluster through the i ntrac luster medium (e.g.. iGunn &: Gott 



19721 : iNulsen et al.lll982l : iFuiita et al.lll999l : lAbadi et al.lll999l ): and/or (2) strangulation, which 



Larson et al. 


1980: 


Bekki et al. 


2002) 



for the MGs in the fields, the mass of the MGs is anti-correlated with the SSFR. This might be 
due to the fact that the more massive MGs have undergone more merging processes than the less 
massive ones and have exhausted more gas and dust in the early precesses. 



7. CONCLUSION 

We assemble a sample of more than 400 merging systems to study their features of star forma- 
tion efficiency. Our sample is drawn from a catalog of ~ 15000 morphologically identified merging 
galaxies. We identify their IR counterparts in the SWIRE data and their optical counterparts in the 
SDSS database. Statistical analyses are performed to explore the SSFR distribution, the relations 
between the SSFR and the pair separation and magnitude difference, and the environmental effects 
on the SSFR. We summarize our main results in the following: 

1. The SSFR distribution of the merging galaxies follows a bimodal distribution. The SSFRs 
of the galaxy pairs that remain unaffected by tidal interactions have a normal distribution; on the 
other hand, the SSFRs of the pairs with enhanced star formation activity approximately follow a 
lognormal distribution. In other words, although most merging systems do show enhanced star 
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formation activity, merging systems with extremely high star formation rates are rare. 

2. There is an anti-correlation between the SSFR and the projected separation of the merging 
galaxies. The SSFR of the pairs also declines when the magnitude difference between two member 
galaxies increases. The pairs with comparable luminosities produce stronger star formation activity 
in total than those with large luminosity contrast; however, in a minor merging pair, the fainter 
component show significantly higher SSFR than the brighter one. 

3. We find a higher ratio of dry to wet mergers in galaxy clusters than that in the field. Dry 
mergers show a higher tendency to populate the clusters than other types of mergers. In addition, 
the SSFRs of wet major mergers are also inhibited in the clusters. 
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Table 1. SWIRE Fields Overlapped with RCS2 



Field 


Center Coordinate (J2000) 


Area (deg^) 




R.A. 


Decl. 




ELAIS Nl 




+55^00"^00'' 


9.00 


ELAIS N2 


16^36'"48" 


+4ldoim45S 


4.45 


Lockman Hole 


10'^45"00'' 




14.32 



Table 2. Ljr/Lz of various types of galaxies based on models in lPolletta et al.l (j2007l ) 



Galaxy type 


Lir/ Lz 


E 


0.11-0.13 


SO 


0.17 


Sa 


0.25 


Sb 


0.42 


Sc 


0.56 


Sd 


1.03 


Sdm 


1.08 


NGC6090 (SB) 


17.25 


M82 (SB) 


20.62 


Arp220 (ULIRG) 


29.57 



-18- 



Table 3. Parameters of the fitting results 



p.d.f. 




a 


normal 


0.65 


0.34 


lognormal 


0.70 


0.43 


normal 


0.65 


0.34 


modified lognormal 


-7.50 


1.14 



Note. — Here ^ and a are the fit- 
ting parameters for the SSFR distri- 
bution of galaxy pairs in the forms of 
the normal and lognormal distribu- 
tion functions. For the normal func- 
tion, ji and a are the mean and stan- 
dard deviation of the SSFR. For the 
lognormal function, /i and a are the 
mean and standard deviation of the 
natural logarithm of the SSFR. 



Table 4. Average SSFRs of bright and faint galaxies in minor mergers and all galaxies in major 

mergers 



CP 


member 


mean SSFR 


Az < 0.05 


minor-faint 


11.7 ± 2.1 




minor-bright 


2.0 ± 0.1 




major 


4.4 ± 0.1 


Az > 0.05 


minor-faint 


11.4 ± 0.8 




minor- bright 


1.9 ± 0.1 




major 


3.8 ± 0.1 
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Table 5. Result of Spearman rank correlation test on Ljr/Lz and Ar^ 



CP 




Ps 


N 


< 0.05 


-0.33 


0.98 


49 


Az > 0.05 


-0.17 


0.95 


132 


0.05 <Az< 0.10 


-0.14 


0.50 


24 


0.10 < Az < 0.20 


-0.04 


0.22 


57 


Az > 0.20 


-0.24 


0.98 


51 



Note. — The parameters and Pg 
are the correlation coefficient and the 
probability of the correlation for Spear- 
man rank correlation test. The parame- 
ter N represents the number of CPs in- 
side each bin of Az. 



Table 6. Comparison between the mergers in the cluster and in the field 





Environment 


Wet merger 


Dry merger 


Mix merger 


Dry-to-wet ratio 




cluster 


4 


1 


5 


0.25 


class 1 


field 


143 


4 


21 


0.03 




fraction in cluster 


0.03±0.01 


0.29±0.16 


0.21±0.08 






cluster 


4 


3 


3 


0.75 


class 2 


field 


139 


12 


17 


0.09 




fraction in cluster 


0.03±0.02 


0.24±0.10 


0.18±0.08 





Note. — In class 1, the galaxies with the Ljr/Lz upper limit lower and higher than 0.5 
are classified into red galaxies and blue galaxies respectively. In class 2, all galaxies without 
24 /xm detection are classified into red galaxies. 
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Table 7. Comparison for wet major mergers in clusters and in field 





number 


mean SSFR 


mean separation (kpc) 


cluster 


4 


2.28±0.42 


13.29±0.74 


field 


109 


1.36±0.04 


13.14±0.07 



6202 (2) 6219 (2) 6223 (1) 6266 (3) 




9031 (2) 9111 (2) 9180 (1) 9289 (3) 




9340 (3) 9352 (1) 9671 (1) 9950 (3) 



Fig. 1. — Examples of merging systems without contamination of stars. The number below each 
image is its original ID in the merging galaxy catalog; the number in parenthesis represents its 
classification in this study: 1 for MG, 2 for CP, 3 for CM. 
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Fig. 2. — Redshift distribution of merging galaxies. 
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Fig. 3.— 



Distribution of redshift difference for galaxy pairs. 
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9136 (0.30) 9160 (0.23) 9339 (0.15) 9351 (0.16) 




9625 (0.38) 9724 (0. 12) 9857 (0.40) 9886 (0. 13) 



Fig. 4. — Examples of CPs with Az > 0.05. The number below each image is its original ID in 
the merging galaxy catalog; the number in parenthesis represents the redshift difference between 
its member galaxies. 
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Fig. 5. — Plots of the 3.6 fim versus the z-band luminosity of galaxies. 
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Fig. 6. — SSFR distribution of galaxy pairs. Middle panel shows a combined model of normal and 
lognormal fits; and the bottom panel shows a combination of normal and modified-lognormal fits. 
In the top panel, dotted and dash dot lines represent the CPs with and without detections at 24 /xm 
respectively, and the solid line represents the overall sample. 
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Fig. 7. — Ljr/Lz as a function of projected separation for CPs with Az < 0.05 (solid circle), 
Az > 0.05 (plus), and MGs (triangle). 
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Fig. 8.— 



Distribution of magnitude difference for CPs. 
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Fig. 9. — Ljji/Lz versus Ar^ for CPs with Am^ < 1.5 (left) and for CPs with Am^ > 1.5 (right). 
Solid circles are the CPs with Az < 0.05; pluses are the CPs with Az > 0.05. 
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Fig. 10. — Ljji/Lz as a function of Am^ for the CPs with Az < 0.05 (sohd circle) and with Az > 
0.05 (plus). The dotted line denotes a dividing criterion between the major and the minor mergers. 
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Fig. 11. — Lir/Lz as a function of Am^; for the minor mergers with Az < 0.05 (left) and those 
with Az > 0.05 (right). Sohd circles are the faint member galaxies; open circles are the bright 
member galaxies. 
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Fig. 12. — Lir/Lz as a function of for the MGs in the field. 



